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Abstract: We consider five-dimensional gravity coupled to a negative cosmological constant 
and a single U{1) gauge field, including a general set of four-derivative interactions. In this 
framework, we construct charged planar AdS black hole solutions perturbatively and consider 
the thermal and hydrodynamic properties of the plasma in the dual CFT. In particular, we 
calculate the ratio of shear viscosity to entropy density and argue that the violation of the 
KSS bound is enhanced in the presence of a chemical potential. We also compute the electrical 
conductivity and comment on various conjectured bounds related to this coefficient. 
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1. Introduction 

Recent years have seen a remarkable confluence of string theory and nuclear physics, with 
efforts to gain new theoretical insights into the strongly coupled quark-gluon plasma (sQGP) 
using holographic techniques [1]. The AdS/CFT correspondence [2, 3] has proven to be a 
powerful tool to investigate the thermal and hydrodynamic properties for certain strongly 
coupled gauge theories [4]. Of course, the gauge theories which are amenable to such holo- 
graphic study are somewhat exotic compared to QCD but one may suppose that certain 
properties of the corresponding plasmas may be universal. The latter suggestion was re- 
inforced by the observation that the ratio of shear viscosity to entropy density seemed to 
be a universal property of the holographic theories yielding rj/s = 1/Att [5, 6]. Further, 
experimental data indicates that this ratio is also unusually small for the sQGP and even 
appears to yield roughly rj/s ~ l/47r [7]. It is now well understood that the holographic 
result rj/s = I/Att emerges for gauge theories described by Einstein gravity as the gravita- 
tional dual. Still this encompasses a remarkably wide class of theories and situations, e.g., 
with various gauge groups and matter content, with or without chemical potentials, with 
non-commutative spatial directions or in external background fields [5, 6]. It is also well 
understood that higher curvature corrections in the gravitational dual will modify this ratio 
[8]-[14]. In fact, it was shown that for certain theories these corrections produce even lower 
values [12, 13, 14] thus disproving^ a longstanding conjecture that rj/s = 1/Att represented 



See also [15]. 
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a strict lower bound for the viscosity of any physical system, i.e., the KSS bound [16]. Still 
one may interpret these new holographic calculations with higher curvature interactions as 
broadening the universality class of conformal gauge theories under study [11, 14]. 

The focus of the present paper is to use the AdS/CFT correspondence to investigate 
how a nonvanishing chemical potential fi effects the hydrodynamics of strongly coupled gauge 
theories. In particular, we consider how rj/s is modified at finite /i. As noted above, with an 
Einstein gravity dual, the result remains ry/s = l/47r, even though individually the viscosity 
and entropy density have a complicated dependence on fi [17, 18]. Hence, fx can only modify 
this ratio through the correction terms appearing from higher derivative interactions in the 
gravitational dual, as we will explicitly illustrate. Since the effects of the chemical potential 
are tied to the higher derivative interactions, it is interesting to examine violations of the 
KSS bound in this context. For example, one might find that the chemical potential limits 
any violations and that the bound is restored with sufficiently large i.e., rj/s > I /Air for 
fi > jJLc- However, we identify a broad class of theories where in fact the opposite result is 
found, i.e., increasing /x only enhances the violation of the KSS bound. 

In principle, studying the effect of the chemical potential on hydrodynamic properties is 
also of phenomenological interest. The higher derivative modifications arc associated with 
corrections emerging from finite Nc and A in the QCD plasma and these may be significant 
for the sQGP [10, 11, 14]. So it is again of interest to determine whether finite ^ enhances or 
suppresses these effects. Unfortunately the relevant chemcial potential for baryon number is 
not expected to be large, i.e., ~ SOMeV or i-i-a/T ^ 0.15 for recent experiments at RHIC 
[19] and so any effects will be limited. However, they may still play a role as the determination 
oi Tj/s becomes more precise in the coming years. 

Turning to the holographic hydrodynamics described by the charged black holes more 
broadly, we also investigate the conductivity, a. It was suggested that the ratio of the conduc- 
tivity to the shear viscosity could obey a bound similar to 77/s [20]. The heuristic reasoning 
behind this conjecture was as follows [20]: in any four-dimensional CFT, we expect rj ~ cT^ 
while a ~ kT where c and k arc basically central charges of the CFT. The first of these is 
related to the total number of degrees of freedom while k is related to the charge degrees 
of freedom. Thus it is natural to expect an upper bound on (Tr^/(rye^) oc k/c, which in 
turn may be related to the weak gravity conjecture of [21]. While this ratio depends on the 
relative normalization of the current and the stress tensor in the CFT, it was also suggested 
in [20] that this relative normalization would not appear in the ratio of the conductivity to 
the susceptibility and so it may be more natural for this ratio to obey a universal bound. 
We extend this discussion to a framework of general four-derivative interactions, as described 
below. 

An overview of the paper is as follows: In section 2, we present our action for five- 
dimensional gravity coupled to a negative cosmological constant and a single C/(l) gauge 
field, including a general set of four-derivative interactions. Further we examine how the 
higher-derivative terms modify charged planar AdS black holes, both the solution and their 
thermodynamic properties. In section 3, we investigate the hydrodynamic properties of these 
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black holes with the four-derivative corrections. In particular, we calculate both the viscosity 
and the conductivity of the dual CFT plasma. Finally, a concluding discussion is presented 
in section 4. We also show in appendix A that one can use field redefinitions to reduce the 
most general four-derivative action to include only the five interactions explicitly studied in 
main text. 



2. Charged black holes in higher derivative gravity 

We begin with five-dimensional gravity coupled to a negative cosmological constant and a 
^(1) gauge field in the following action: 



I = 

21 



p I d'xV^ g + _ + '^e^'^'^'AaFbcF^ + [cRabcdR''"^ (2.1) 
+C,RabcdF'''F'^ + CsiF^f + c,F^+ c, s'^'^'^'AaRbcfgRdJ')' 



where F'^ = FabF"-^ and = F^'^F'^cF'^aF'^a- As well as the conventional Einstein and 
Maxwell terms, our two-derivative action also includes the Chern-Simons term proportional 
to e^^'cde^ which naturally arises in five-dimensional supergravity [22]. The above action (2.1) 
also contains a general set of four-derivative interactions. We will treat these terms in a 
perturbative framework where each of the the dimensionless coefficients C; <C 1. As discussed 
in [14], it is natural to expect that each of these coefficients is suppressed by a factor of 
which we are assuming is very small. We demonstrate in appendix A that within 
this perturbative framework, one can use field redefinitions to reduce the most general four- 
derivative action to include only the five interactions appearing above in (2.1). 
The metric equation of motion arising from (2.1) is 

1116 

Rab-^Rgab = ^FacFb'--F^gab + J^gab (2-2) 
+L'^Ci_ Rede fR"^^^ gab - '^R{a\edeR\bf'^'^ + 4V''V'^i?c(a6)d^ 
+L^C^ (^^RcdefF'^'F^^gab + 3R'''\aFb)eFed + 2^=^^ (F,(„F,,)d) j 

+21^0,8""'^^ {R'\b)e}^gF^ + 2F,dVeRflb)) 

while the vector equation of motion is given by 

VbF^" + Ke'''"^^FbcFde = -AL'c^Vb{R''^"^F^d) (2.3) 
+SL\,Vb (f^F'"^) + SL^c^V^ (F'\F\F'^d) - L^c,e''''"''RbcfgRdJ' ■ 



As discussed above, we will solve these equations of motion perturbatively in the coefficients 
C; of the interactions in the four-derivative action. 
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2.1 AdS/CFT dictionary 

With a negative cosmological constant, the gravitational theory described by (2.1) naturally 
has an AdS5 vacuum and is dual to a four-dimensional CFT. In this holographic context, the 
bulk vector field will be dual to the current generating a global U{1) symmetry in the CFT. 
In all, the action (2.1) is characterized by seven dimensionless parameters: L^/£p^, k and 
the five coefiicients Cj. The AdS/CFT correspondence then relates each of these gravitational 
couplings to various parameters that characterize the dual field theory. For example, the 
holographic framework relates the two central charges, a and c, of four-dimensional CFT to 
[23, 24] 

cf^3c — a\ Ic — a , 
1-- , Ci ~ . 2.4) 

As described in [14], a key assumption in working with the effective action (2.1) is that 
the five-dimensional gravity theory is described by a sensible derivative expansion. That is, 
we are implicitly assuming that couplings of the four- and higher-derivative interactions are 
systematically suppressed by powers of the Planck length over the (bare) AdS scale, ip/L. 
In particular then, we expect that Ci oc (.p^ /L? <^ 1. From the perspective of the AdS/CFT 
correspondence then, we are restricted by (2.4) to consider CFT's for which 

c ~ a > 1 and |c - o|/c < 1 . (2.5) 

Further, our assumption about the derivative expansion in the gravity action then restricts 
the size of the field theory parameters related to the four-derivative couplings, i.e., the CFT's 
of interest should have the corresponding parameters being suppressed by inverse powers of 
the central charge c. 

Turning to the other couplings, it is natural to consider k and C5 together since they 
both appear in interactions proportional to e"^'^'^. Both of these Chern-Simons-like terms are 

not invariant under 'large' gauge transformations and as a result, in the present holographic 
context, they play the distinguished role of determining anomalies for the global U{1) sym- 
metry in the dual CFT [25, 26]. In our present analysis, we leave these coefficients to be 
arbitrary constants but we should also note that their precise values are irrelevant here since 
these terms play no role below in determining the geometry or thermal properties of our 
background solutions. 

In contrast, the interactions parameterized by Ca, C3 and C4 all play a role in our perturba- 
tive analysis of the charged black holes, as well as Ci. From the point of view of the dual CFT, 
C2 characterizes the form of the three-point function of two currents with the stress tensor 
[27]. Similarly, C3 and C4 provide two independent couplings in the four-point function of four 
currents. Again, we leave all of these coefficients arbitrary in our general analysis. However, 
in the context of a supersymmetric theory, a special case arises for the i?-symmetry current 
which is in the same supermultiplet as the stress tensor. In this case, all of the corresponding 
CFT couplings will be proportional to the difference of the central charges [27], as found for Ci 
in (2.4). The holographic dual of such a supersymetric CFT is an A'^ = 2 supergravity theory. 
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While the latter may be gauged or ungauged depending on the details of the CFT, the dual 
of the i?-symmetry current is the particular U{1) vector appearing in the five-dimensional 
graviton supermultiplet. In this framework, supersymmetry dictates the form of the four- 
derivative corrections to the leading supergravity action and so all of the relevant couplings 
Ci are again related [26, 28]. This supersymmetric setting will be of particular interest in our 
discussion in section 4. 

We should note that there is one other (dimensionless) parameter implicit in our analysis, 
which can be described as the relative normalization of the gauge and gravity kinetic terms 
or alternatively as the ratio of the five-dimensional gauge coupling to, say, the AdS scale. 
One can see there is an issue here since in the conventions used in (2.1), the gauge field is 
dimensionless while a conventional gauge connection should have units of energy or inverse 
length. Hence, we should scale the gauge field by some appropriate scale, = Aj^. With 
this choice, the Maxwell term in (2.1) becomes J d^x^/^F^ where the five-dimensional 
gauge coupling is given by = 2£p'^/Ll. In any particular setting, one is typically guided 
by the details of the AdS/CFT correspondence or the string theory construction to give the 
proper normalization of the gauge field, i.e., choosing the scale — for example, sec [17, 29]. 
For simplicity, in the following we make a particular convenient choice for L^, but of course 
it is a straightforward exercise to reinstate a general in our results. 

To close this subsection, we observe that typically in supergravity actions, the gauge 
kinetic terms will couple to various scalars. Prom the dual CFT perspective, such coupling 
would indicate a nontrivial three-point function mixing two currents with some scalar oper- 
ator. Hence, from this point of view, the action (2.1) is not the most general since we are 
making a special choice for the form of the vector kinetic term in the two-derivative action. 
Beyond this choice, we note that while we are also dropping any possible scalar couplings in 
the four-derivative interactions, such couplings would only contribute at the next order in our 
perturbative expansion [14] . 

2.2 Charged black hole solutions 

We consider charged planar black hole solutions with the following ansatz:^ 

^^2 ^ _!:!ffirfi2 + J^dr' + ll^dx' + dy^ + dz') , (2.6) 
j^z r-^gyr) 

At = h{r) . 

The leading order solution (of the two-derivative equations of motion) may be determined to 
be: 

/. = «.= (l-|)(l + ^-|i). (2.T) 

= 5<?i' (4 - ^) = ^ ■ (2.8) 

^Charged black hole solutions with spherical horizons were constructed for a general four-derivative action 
in [30]. 
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Here Tq denotes the position of the (outer) event horizon. There is also an inner horizon at 




(2.9) 



The solution is characterized by the charge density, which is given by {*F)^y.^ = Q, and the 
mass density, which is proportional to M = VQ+q^/rQ. Implicity, we have fixed the integration 
constant in ho such that the gauge field vanishes at the horizon, as required by regularity.^ 
This leading order solution is extremal with q^/vQ = 2 for which (2.9) shows the two horizons 
coincide, i.e., = r„. With q'^/rf^ > 2, > rf^ and the solutions actually describe the 
same set of noncxtremal black holes as with q^ /rf^ < 2 but with r„ and r_ exchanging roles. ^ 
Solutions where ratio of charge to mass densities exceeds that in the extremal black hole 
{i.e., Q^/M^/^ > 8/\/3L^) are found by allowing to become negative but, of course, such 
solutions all contain a naked singularity at r = 0. 

Now we wish to construct pcrturbative solutions to first order in the C;. We maintain the 
ansatz (2.6) and parameterize the perturbative solution as 



where F{r), G{r) and H{r) are 0(ci) corrections. It is then straightforward to solve the 



^To see this, note that the event horizon in (2.6) is the Killing horizon where |9t|^ = 0. However, as a 
Killing horizon, it also contains the bifurcation surface which is a fixed point of the Killing fiow, i.e., dt = 
on the bifurcation surface, as opposed to the previous null condition [31]. Hence if the gauge field A is to bo a 
well defined one-form, then At must vanish there. This is, of course, the Lorentzian analog of the topological 
constraint which arises more intuitively for the corresponding Euclidean black hole. 

^This symmetry between ra and r_ is readily seen by noting that (2.9) comes from demanding the vanishing 
of the second factor in /o, i.e., Tq ri + r^r^ — = 0. 



f{r) 
9{r) 
h{r) 



/o(r)(l + F(r)), 
/o(r)(l + F(r) + G(r)), 
ho{r)+H{r), 



(2.10) 
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equations of motion (2.2) and (2.3) to first order:^ 

4 2 8 / 2 \ 2 

6/oF(r) = 2(2ci-35i) + ^/i + 12%(^2-26ci-12c2) + 12% 1 + %) c, (2.11) 

r r ^ V ^0 / 

2 4/ 2 \ 4 

+8^ ( 1 + % ) (5ci + 6C2) + ^2 (17ci - 24(c2 + 6C3 + 3C4)) , 
G'(r) =5i-^^(13ci + 12c2), (2.12) 



H{r) = h,- Vs^^h, - 8^3^ (1 + ^ 



, 6 ' '^^ 

+;^-^(48c2 + 144c3 + 72c4 - 13ci) , (2.13) 

where /i, 51, hi and /i2 are (dimensionless) integration constants. 
We fix the integration constants as follows: 

• The background metric for the dual CFT can be extracted from the asymptotic be- 
haviour in the black hole metric (2.6) as 

dslp^ = -foodt'^ + dx'^ + dy"^ + dz'^ (2.14) 

where we defined foo = f{T~ ^ 00). Hence to fix the speed of light to be one in the dual 
gauge theory, we require that foo = 1. From (2.11), we find 

91 = \c, . (2.15) 

Note that this fixes the asymptotic behaviour g{r — 00) — 1 + 2ci/3 which reflects the 
fact that, as noted in [14], the AdS scale of the background geometry is perturbed to 
be 

^ = ;^(l + |c.). (2.16) 

when Ci is nonvanishing. 

• For simplicity, we require a regular F{r) and fix the position of the event horizon to 
remain at r = Tq. That is, we require that foF{r = To) = 0. Again from (2.11), this 
fixes /i to be 

h = -2 (8ci - 351) + 4\ (62ci + 24c2 - 3/12) (2.17) 

2\ 2 

% ) (69ci + 24(c2 - 6C3 - 3C4)) . 



■'Our approach was as follows: Examining the linear combination of the metric equations (2.2) proportional 
to G't — f /g Cr (where Ci, is the Einstein tensor), one finds a first-order linear ODE for G(r) which is readily 
soluble. Given G(r), the t component of the vector equations (2.3) is easily solved for H{r). Finally with the 
solution for those two perturbations, F{r) can be determined by solving the first-order linear ODE coming 
from the r equation. 
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• We require that At vanishes on the horizon — as described in footnote 3. Setting 
H{r = To) = in (2.13), we fix to be 



hi = V3^{h2 + 8c,) + ^^(13ci - 24(c2 + Gcg + 3c^)) . (2.18) 

Finally we may use the remaining freedom to require that the charge density is fixed 
as in the leading order solution, i.e., {*F)^y^ = Q. The perturbed vector equation of 
motion (2.3) can be written in the form V^X''" = for the appropriate antisymmetric 
tensor Xab- In the perturbed solution, i*X)^y^ is a constant independent of radius and 
it is natural to define the charge density to be {*X)^y^ = Q. This allows us to fix /12, 
which is most simply done by examining this constraint for asymptotic r where 



lim {*X)^y^ = lim 



r^3 



r— >oo 



~ (^l + ^gi + h2 + 8c,^Q. (2.19) 

Hence we fix 

h2 = -^91- 8c, = -^ci - 8c, , (2.20) 
where in the last expression, we substituted for gi as in (2.15). 

2.3 Black hole thermodynamics 

For the thermodynamics of the above charged black holes, let us begin by first reviewing the 
results for the leading order solution, (2.7) and (2.8).® The temperature of the dual CFT is 
precisely the Hawking temperature calculated as the inverse of the periodicity of time in the 
corresponding Euclidean solution: 

'l-:&V (2-21) 



7rL2 V 2r6 



Note that the temperature vanishes for the extremal black hole with jr^ = 2. 

Next we would like to consider the chemical potential of the system which is related to 
the asymptotic value of the potential At- However, as discussed in section 2.1, we must scale 
the gauge field by some appropriate scale, A^ = L^, A^, to produce a chemical potential with 
the appropriate units of energy. The chemical potential then becomes 



®Thc thermodynamics of charged AdS black holes has been well studied [32, 33], of course, but the focus 
was on solutions with spherical horizons. 
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In any particular setting, one would be guided by the details of the AdS/CFT construction 
or the string theory construction to give the proper normalization of the chemical potential. 
However, for simplicity in our general analysis, we will make the convenient choice 

L^=ttL (2.23) 

in the following. Of course, it is a straightforward exercise to reinstate a general L* in the 
following calculations. Note that with the preceding choice, we may write 

'T . . ■ (2.24) 

Further in the extremal limit T = 0, the horizon radius remains finite with ro = 7rL^/i/\/6. 
It will also be convenient to denote the ratio of the chemical potential to the temperature as 

This formula may be inverted to yield 





where the last expression is a Taylor series for small fi7 

To proceed further, we apply the standard path integral techniques [34] in which we 
identify the Euclidean action Ie = W/T, where W{T,iJ,) is the Gibbs free energy, i.e., the 
thermodynamic potential in the grand canonical ensemble.^ To calculate the Euclidean action, 
as well as the bulk action (2.1), one includes the Gibbons-Hawking boundary term [34] and the 
appropriate boundary 'counter-term' action [23, 35]. Alternatively, one can use background 
subtraction and consider the difference in the (bulk) action for the charged black hole and 
for AdS5 with a constant gauge potential. Further since we are considering planar black 
holes (2.6), the spatial volume in the dual CFT is infinite and so we divide our result for any 
extensive quantities by a regulator volume Vx and work with the corresponding density. We 
do not go into the details of the calculations here but only present the final result for the free 
energy density: 

+ ^f^'r! ] . (2.27) 



^We provide such a Taylor series expansion for all of our results with an eye towards the fact that /xb /T is 
small at RHIC. 

®Onc could also consider the microcanonical ensemble with a fixed charge density by making the standard 
Legendre transform to the Helmholtz free energy: F{T, riq) = W{T, n) + J cPx Uq jjb. In the AdSs language, this 
corresponds to adding an additional boundary term to the action which ensures that the appropriate boundary 
condition corresponds to fixing the (radial) electric field, rather than the gauge potential, at asymptotic infinity 
— for example, see [29, 32]. 
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While in principle we could use (2.24) to express the free energy density entirely in terms of 
T and /i, the last expression above with w {ro{T, n), fi) is sufficient for most calculations. In 
particular, the standard thermodynamic identities yield the entropy density and the charge 
density: 



s = — 



rin 



dw 




2'Krl 


&T 






dw 






djjL 


T 





(2.28) 
(2.29) 



Note that the result for the entropy density matches the expected result for the Bekenstein- 
Hawking entropy of the black hole horizon. Using the previous expressions, we can also 
express these quantities in terms of the temperature and chemical potential: 



s = 



3"^ 




(2.30) 




(2.31) 



where in both cases, we have also given a Taylor series for small jl. Hence for large T (or small 
n/T) our holographic model yields s ocT^ and Uq oc /xT^, as may have been anticipated. We 
can also combine the above expressions, (2.28) and (2.29), to calculate the energy density: 



PE = W + Ts + IXUq 



3 rt 



2 Ip^L^ 



1 + '- 



3 M 



2 e/L^ ' 



(2.32) 



We now turn to the first order solution (2.10) which takes into account 0(ci) terms in 
the equations of motion. The temperature of the perturbed black hole becomes 



^a-|g)(i + ^('-c) + iGM 



7rL2 



5 / 31 

'-3^^-2;^^ + y^^ + '^^^ 



(2.33) 



,2\ 2 



^ (9ci-4c2-24(2c3 + C4)) 



Note that the corrections shift the condition for the extremal limit T = to be 

^2 



Q 



2 [1 - 48 (ci - 2(2c3 + C4))] . 



(2.34) 
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The asymptotic value of At determines the chemical potential, as in (2.22). The modified 
result is given by 



V3( 



(2.35) 



1 

I - -Ci + 8c2 + -g 
L ^ ^0 



13 



Ci - 8c2 - 24 (2c3 + C4,] 



where we have chosen L^, as in (2.23) and hi is fixed as in (2.18). 

The free energy density is most easily calculated using background subtraction, as de- 
scribed above, with the final result: 



w = — 



3-^ 

2ip 



19 



113 



1 + + ^ 1 - — ci - 32c2 + ^ hr^i + - 12(2c3 + C4 



23 



1 + 13ci + (1 + yCi)A^ 



26 



(2.36) 



+ ( 1 + y ci + 24(c2 + 2c3 + 



6 



108 



(1 - 15ci)/i6 + 



Now using the same thermodynamic identities as above, we arrive at the following expressions 
for the entropy and charge densities: 



27rrg 

4:ip 



l + 8ci -4(7ci + 6c2)\ 



I + W1 + -/22 +4 I + ./1+ ^: 



78-2^2 + 6-^^ + ^^ 



'l + lA^ 



1 + 13ci + 



11 A At . 



n„ 



l + ^/l + -/l2 



+ 



9 33 - An 
33 + 46/i2 + ^ 



2r,2 



£ ^ 



-, 11 

1 H Ci + — 

3 3 



2 ^ 26 

1 + — Ci + 24(c2 + 2c3 + C4) 




. 73 , 
l-yc,) + 



(2.38) 



Note that the first expression for the entropy density matches precisely the result found using 
Wald's formula for higher curvature theories [36]. It is interesting to observe that when the 
entropy is expressed in terms of T and fi, it becomes independent of C2, which appears in the 
original 'geometric' expression for the entropy. In contrast, all of the coefficients, c^, Ca, C3 
and C4, appear in the charge density. We also note here that just as for the leading order 
results, one finds pE = —Sw when the first order corrections are included. 
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3. Holographic hydrodynamics 



We now turn to computing the shear viscosity and the conductivity of the holographic plasma 
represented by the charged black holes in the previous section. We follow the approach 
of expressing the transport coefficients in terms of field theory correlators using the Kubo 
formula and then calculating these correlators with holographic techniques [37]. However, 
as we are working with the higher curvature interactions in (2.1), we must generalize these 
standard calculations. The analogous calculations for a particular four-curvature interaction 
first appeared in [8] and the latter are readily adapted to the present case. Our presentation 
also builds on the recent work of [6] which focussed on the hydrodynamic limit of low frequency 
and momenta, showing that there is a simple relation between quantities computed in the 
membrane paradigm approach and those calculated using AdS/CFT. In particular, the shear 
viscosity of a field theory with a gravity dual may be related in a simple fashion to the 
membrane coupling constant of a certain minimally coupled scalar in the dual gravitational 
background. Our calculation closely parallels the discussion of [38] which presented a general 
framework to calculate the shear viscosity for higher curvature theories. 

Before proceeding, we make a change of coordinates u = r^jr^ which is more readily 
adapted to the hydrodynamic calculations. With this coordinate choice, the background 
solution (2.6) becomes 

= + -^Au^ + 4-(rfx2 + V + dz^) , (3.1) 
At = h{u) , 

where the leading order solution (2.7) and (2.8) now takes the form: 

fo = 9o = {l-u)(l + u- '^uA , (3.2) 



^0 = (1 - n) . (3.3) 

Of course, one must also make the appropriate substitution in the perturbative solution given 
by (2.10)-(2.13). A simplifying feature, however, is that with our choice of the integration 
constants described above, in particular for /i, the event horizon remains fixed at r = ro or 
u = uq = 1 m. the perturbative solution. Of course, in terms of the new radial coordinate, 
the asymptotic boundary corresponds to n = 0. 

3.1 Corrections to r]/s 

Kubo's formula relates the shear viscosity to the low frequency and zero momentum limit of 
the retarded Green's function of the stress tensor in the CFT 

G%,xy{^,'^ = 0) = dtd^^"^e{t) {[T,y{x),T,y{m . (3.4) 
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Concretely one has 

rj = - lim - Im {uj, k = 0) . (3.5) 

uj^O CO 

The retarded Green's function may be computed using the prescription first set out in [37]. 
Translating the calculation of the correlator to a holographic one, one first finds the effective 
action for the metric perturbation hx^{t,u) = J (/>fc(it)e~*'^*+''^^. Evaluating the action 
(2.1) to quadratic order in the fluctuations ^fc(«) yields 



j{2) 



1 f d'^k 

+D{u)<t>kcf>-k + E{u)4cf>'Lk + F<t>l<t>U) + ^- (3-6) 



This form of the effective action originally appeared in [8], where the effect of certain 
terms were considered, but this general form will arise for any action involving any powers 
of the curvature tensor (but not derivatives of the curvature^) and so appears again in the 
present context with the action (2.1). If we consider an action where the higher derivative 
terms come coupled through some parameter 7, then the two functions E and F are 0(7). 
Following [8], we have also added a generalized Gibbons-Hawking boundary term /C in (3.6), 

/C 



with 



^^/(l^^^^+^^ + ^^)l-- ^''^ 

K, = -Acl>',cf>^k K2 = ~4>'k^U (3-8) 

K3 = E{pi^'k + 2po(l>k) (t>'-k 

The first term Ki is essentially the contribution of original Gibbons-Hawking term while 
K-2 and are new 0(7) contributions. In the term K^, the coefficient functions PQ-,pi are 
defined in terms of the linearized equation of motion for (f): 

(t>"+pi(t>' + poct) = 0{j). (3.9) 

With this boundary term (3.7), the variational principle is valid up to 0(7^) [8]. 

Given the effective action, let us proceed in trying to compute the shear viscosity. First, 
we integrate by parts to rewrite the action in a more symmetric fashion: 

(2) _ 1 f d'^k 



21 



1 r d^k 

p J ^,du {{B-A- F' /2)4>'^4>'-k + Ect>l4>'Lj, 

_ (c - A')72)</'fe<^-fe) + it . (3.10) 



The integration by parts eliminates Ki and K2 in the boundary term and /C is given by: 

u=l 

(3.11) 



u=0 



generalization to include derivatives of the curvature appears in [38]. 
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At this point, it is convenient to define the (radial) canonical momentum for our effective 
scalar as: 

'^^^ 1 

n^^^) ^j;^ = 7^(^^-^- ^V2)0U«) - (E^M)') , (3.12) 

where in the variation we regard 0" as {(/)')'■ The scalar equation of motion then takes the 
simple form 

duIlkiu)=M{u)Mu), M{u) = -^{D-{C-Ay/2). (3.13) 

together with the definition of the canonical momentum (3.12). 

To compute the retarded Green's function, we now evaluate the effective action on-shell, 
which reduces to a boundary term using the linearized equation of motion (3.13) 

-^^i::;- (3-14) 



(27r) 



The retarded Green's function is then given by the flux factor [37] evaluated at the asymptotic 
boundary 

G^^ _(a;, k) = - lim ^ , ^^f"^ , , (3.15) 

^y^^y^ (/)fe(n)0_fe(n) 

where the factors (f)k{u)(f)-k{u) ensure the appropriate normalization for the Green's function. 
Implicitly, this expression is also evaluated on <p{u) with inf ailing boundary conditions at the 
horizon. In the present case, the flux factor reduces to 

2J^k = nfc(/)_fc + (C - A')(t>k<l)-k + E(^'i:<t>'_k + 

= Uk^P-k + {C- A!)(t>kct>-k + Epo<p',,<j)_k, (3.16) 

where on the second line we have used the lowest order equation of motion (3.9) for (pj.. We 
see that the flux is given almost entirely by the canonical momentum term. However, at this 
point we note that, since according to (3.5) the shear viscosity is given by the imaginary part 
of the Green's function, then the second term will not contribute as (pk'P-k is real. It turns 
out that we can discard the third term as well since it is of 0(a;^), as we now explain. 

It is an important point that in general the effective mass M{u) in (3.13) is 0{u}'^) and 
therefore can be set to zero in the low frequency approximation. Consider setting 4>{u) to a 
constant, in which case the corresponding radial momentum (3.12) automatically vanishes. 
The equation of motion (3.13) must be satisfied, since a constant (p{u) simply corresponds 
to a rotation and rescahng of the x,y coordinates. Therefore M{u) must be 0(w), but time 
reversal invariance demands that it must be proportional to w^. We conclude that on general 
grounds we must have M{u) = 0{lo'^), and therefore can be set to zero in the low frequency 
limit, which is taken in calculating the shear viscosity via (3.5). In particular, with regard 
to the flux in (3.16), the third term is proportional to the mass term of the lowest order 
equation of motion and so is 0{lo^). Hence this term is also irrelevant in calculating rj in 
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the low frequency limit. We conclude that the only relevant piece in the flux is the canonical 
momentum term, and so 

77 = - lim -ImG^„„„(w,k = 0) = lim ^^'^\ ,. (3.17) 

oj-^OLO ^^'^^^ ^ u,oj^O lUJ(j){u) ^ ' 

Here H('u) = n{j^ k=o}(^) ^'^id ^(n) = '?!'{a;,k=o} ("") • Further, in this limit, the equation of 
motion (3.13) for the canonical momentum is simply 

a„Hfc(u)=0, (3.18) 

i.e., nfe('u) is independent of the radius. Therefore we are free to evaluate the value of n(u) 
in (3.17) at any radius and in particular, it can be evaluated at the horizon. 

Hence, the final ingredient to evaluate (3.17) is to determine cotplu). Now to fix the 
fluctuations, we must impose infalling boundary conditions on (pki^) at the horizon u = uq. 
These together with regularity at the horizon imply [6]: 



du(t>{uo,t) = -iu 
dl4){uo,t) = -iuj du 



-9tt 



-gtt 



dl<p{uo,t) = -iudl' 



-gtt 



(t>iuo) + 0(a;') 



Uo 



Following [6], keeping LO(j){u) and Ii{u) constant in the low frequency limit, the definition of 
n implies 

u(j)'{u) = j{Ci{u)(j)"{u) + C2iu)4>"'{u)) (3.19) 

with some functions Ci, C2 whose detailed form is irrelevant to our purposes. Working pertur- 
bativcly in 7,-'^'' one performs a split 0(ii) = (poiu) +7<^i(n). Then, to lowest order in 7, (3.19) 
yields the solution: u}(j)Q{u) = ujcpo^O), i.e., iO(po{u) is also constant in the low frequency limit. 
At the next order in 7, the equation of motion for (jj(j)i{u) then also reduces to u;(l)[{u) = 
and again with the solution: LO(f)i{u) = a;(/>i(0). We conclude that (jj(f){u) = uj(l){0) to leading 
order in the low frequency limit. 
Hence we arrive at the result 

V = lim = (k2(«o) + «4(uo)) . (3.20) 

where in the second expression, we have evaluated the ratio at the horizon u = uq and defined 
the quantities 

/ X / guu{u) f . N , F'{U)\ i TPt \( I gnu{u) 

H2{u) = 4/ -—- A[u) - B[u) H — , K4,{u) = E{u) ' ' 



gtt{u) V 2 y I V V dttiu) 



(3.21) 



^"For certain higher curvature actions, e.g., Gauss-Bonnet gravity, the equations of motion for are still 
second order in derivatives, in which case Ci = C2 = and our conclusion follows immediately. 
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The indices on K; indicate the number of derivatives appearing in the corresponding terms in 
(3.10). Note that up to now our discussion of calculating rj has been completely general and 
this approach applies to any higher derivative action involving powers of the curvature tensor 
(but not derivatives of the curvature). We have verified that this approach reproduces the 
known results in the literature for theories containing four-curvature [8, 9] and two-curvature 
[12, 13] interactions. 

Now we specialize the discussion to considering the action (2.1). Using the action given 
in equation (2.1), the functions A,B,E,F turn out to be: 



= l + 4ci(/^(u)-/o(ix)) 



A 








B 


3. 






E 








F 


V^5"" 



1 



4 /o(u)^-2uV^(t.)^ 
S"' fo{u) 



4ci 

16cJ^(ii) 



where we have used that to lowest order in C;, f{u) = g{u) = fo{u) 
in hand, it is straightforward to obtain the shear viscosity: 



(3.22) 

With these expressions 



^=^(l-8c.(i^'-/^)). 
Combining this with our previous result for the entropy density (2.37), we arrive at 



s 



^ . l-8ci+4(ci+6c2)\ 
47r V rX 



(3.23) 



(3.24) 



This agrees with the well established result when q = [12, 13]. Note that c-^ and C4 do not 
appear in this expression. We can combine the above with (2.26) to express the ratio in terms 
of /2 = n/T: 



1 
s 



1 

Ait 
1 



16|Lt^ (ci -I- 6c2 



3 1 + ^1 + 2/i2/3 



1 - 8ci + 



1 - 8ci + ^(Ci + 6c2)/i2 1 1 - 1/^2 + 0(^4^ 



(3.25) 



We note that the infalling boundary conditions are modified for the extremal black holes 
and so in principle, our computation would need to be modified in this case. Nevertheless 



^^We should point out that to obtain this result the corrections to the lowest order background play no 
role other than defining the temperature and the chemical potential. This could have been anticipated since 
the corrected background only comes in at the two derivative level, but there on general grounds the shear 



viscosity has the universal form rj = V3/{2lp) where V3 = {qx 



N3/2 



[37]. 
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we can consider our result (3.24) in the extremal limit where i.e., q'^/r^ 2 + 0(ci). This 
limit yields at T = 0: 

^ = (1 + A8c,) . (3.26) 
s 47r 

Notice that the leading correction is now independent of Ci unlike the T ^ case. 
3.2 Conductivity and higher derivative terms 

We now turn to computing the DC conductivity in the pcrturbative background corrected 
by the four-derivative interactions. This can be obtained by using a Kubo formula similar to 
the one for the shear viscosity. Let us define 

G^^(a;, k = 0) = -ij dtd^e'^'eit) ([J^(x), J,(0)]) ., (3.27) 

where is the CFT current dual to the bulk gauge field A^. Then the DC conductivity is 
given by: 

(7 = - lim *- Im G^^iu, k = 0) . (3.28) 

Here the factor appears in the prefactor so that a corresponds to the conductivity of 
the current dual to the properly normalized potential Further, in order to interpret the 
result as the 'electrical' conductivity of the plasma, we imagine coupling the CFT current 
to an external or auxiliary vector field, following [20, 39]. This auxiliary vector gauges the 
corresponding global U{1) symmetry in the CFT with a (small) coupling e. Then to leading 
order in e, the effects of the auxiliary vector are negligible and the conductivity can be 
determined from the original CFT alone. The same result can also be related to the thermal 
conductivity [17] which determines the response of the heat flow to temperature gradients, i.e., 
T^i = —K-T diT in the hydrodynamic limit. The full expression for the thermal conductivity 
can be written as [17] 

\ 2 rri 

s M \ 



The computations are most conveniently performed within an effective action approach, 

as in the previous subsection. Since the At component of the bulk vector is nonvanishing 
in the background (2.6), the perturbations Ax can couple to the shear mode graviton, i.e., 
metric perturbations of the form hxi- However, gauge invariance imposes a relation between 
the two sets of perturbations which we use to integrate out the hxi and obtain an action that 
involves only the Ax fluctuation. 



^^We can use the leading order result found from (2.21) here since q only appears in the correction term in 
(3.24). 

^^The full relativistic expression for the heat flow also includes a contribution proportional to the pressure 
gradient [17]. 
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Starting with the action (2.1), we compute the quadratic action for the perturbations 

^/ife(u)e— (3.30) 

To begin, we consider only the leading order background, (2.7) and (2.8), setting c, = 
in the action. We would like to set the perturbation to zero as a gauge choice. The 
corresponding component of Einstein's equations then becomes a constraint which yields: 

9xx t'k = -A[ ak . (3.31) 

Plugging this constraint back into the effective action along with hk{u) = 0, the quadratic 
action takes the simple form: 

= ^ / -^du {N{u)a'ka'_k + M{u)aka^k) , (3.32) 

where 

N{u) = -^Mu), M{u) = - -A'^iuf . (3.33) 

The equation of motion is solved near the horizon with the ansatz 

afc(«) = CUuT (3.34) 

with usual. The infalling boundary condition corresponds to choosing the 

minus sign. The equation of motion for can be re-expressed as 



dujk{u) = j^M(u)ak{u) (3.35) 

where, as in the previous section, we have defined the radial momentum for the effective scalar 

df^'^ 1 

Uu) = ^ = -T^N{u) a',{u) . (3.36) 
The condition of regularity at the horizon u = uq corresponds to setting [6] 



jk{uo) = lim ^^^y^.[^ak{uo) + 0(u;^) , (3.37) 
"-*"o ip-" V gtt 

where we are expanding in small u with the zero-frequency limit of (3.28) in mind. Next one 
evaluates the on-shell action to identify the flux factor, which can be written as simply 

'^^k = jk{u)a_k{u) . (3.38) 
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The Green's function (3.27) is given by evaluating the flux with the appropriate normalization 
at the asymptotic boundary. The DC conductivity (3.28) is then given by a formula analogous 
to (3.17) for the shear viscosity, 

Im[ife(«)a_jt('u)] 



lim Im 

U,U!—*0 



2J^fc 



lim 



(3.39) 



k=0 



_afe(u)a_fe(u)Jk=o 

where it is convenient not to cancel the factors of a_jk(n) in the final expression, as will 
become apparent below. The key difference between the present case and the computation 
of the shear viscosity is that neither nor u is independent of the radial position, even in 
the low frequency limit. As is evident from (3.33), the effective mass M{u) no longer vanishes 
in this limit and so the equation of motion (3.35) still produces a nontrivial flow in the radial 
ridection. However, if we apply (3.35) in examining the radial evolution of the numerator in 
(3.39), we find 

d 



du 



Im[jfe(^i)a-fe(^i)] = Im(/i(M)ofea_fe + f2{u)jkj-k) = 0. 



(3.40) 



Notice that this result does not rely on the fact that we are taking a low frequency limit. 
Therefore we are free to evaluate Im[jfc(n)a_fc(n)] at any radius, e.g., at the horizon. However, 
at the horizon, jk is constrained by the regularity condition (3.37) and so we may write 



4^ '^"^ m) 



k=0 



where we have defined 



4{u) = -N{u) 



-9tt 



and A/fc(u) = ak{u)a-k{u) . 



(3.41) 



(3.42) 



The quantity N{u) is real and so independent of uj up to 0{uj'^). This also means that to this 
order, M{u) is completely regular at the horizon, since the logarithmic divergence of ak{u) 
there is always accompanied by a factor of iu/T. Therefore, in computing M, we are free 
to solve ak{u) imposing regularity at the horizon and setting u to zero, which simplifies the 
calculation considerably. In the leading order background (3.1)-(3.3), the solution is easy to 
obtain: 



ak{u) = ajt(O)- 
The leading order conductivity then follows: 



1 



2^ 



(2 - 3u) 



1 + 



(3.43) 



a = 




(3.44) 
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using our previous formulae (2.23)-(2.26). 

Extending these calculations to work at first order in the four-derivative couplings C; 
is straightforward. One follows the same steps as above. That is, one first computes the 
quadratic effective action for afc,tfc,/ifc and obtains a constraint upon setting = 0. Sub- 
stituting the constraint back into the action and keeping terms linear in C; one still gets an 
action of the form (3.32). In the subsequent steps to calculate a, the equation of motion for 
ak is technically harder to solve and so we only present the results to leading order for small 
p,, i.e., to order 0{p?) or 0{q^). The conductivity then turns out to be 



e L,ro (_^_^ 16c2 + ^(-3 + 68ci + AOc^ + 96(2c3 + c^) 



L 



(3.45) 



o 3^ I 1 + + ^^^2 - - -^(15301 + 112c2 + 192(2c3 + C4))] + 0(A^^ 

Now following [20], we examine the ratio 



ctT^ 4_2 10 

2 = 1 - q/^ - + I6C2 

?7 3 3 



-/i=^[13ci + C2 + 12(2c3 + C4)] + 0(/i^ 



(3.46) 



As described in the introduction, [20] suggested that the simplicity of the leading result at 
p, = 0, i.e., 1, may indicate that this result is universal. Further they argued that this 
leading behaviour may then represent a universal upper bound for this ratio. The above 
result certainly indicates that aT'^/{r]e'^) < 1 even for / 0. However, the leading order 
result is also modified by the four-derivative couplings and so the upper bound conjectured 
here may also be violated depending on the precise values of these couplings, similar to what 
was found for the KSS bound in [12, 13]. 

In [20], it was also noted that the ratio above depends on the relative normalization of 
the current and the stress tensor in the CFT but they further observed that this relative 
normalization does not appear in the ratio of the conductivity to the susceptibility. Hence 
[20] suggested that it may be more natural for the latter ratio to obey a universal bound. 
To examine how the higher derivative couplings effect this ratio, we must first calculate the 
susceptibility S: 



S(T)^ 



dnq 



(3.47) 



where Uq is defined in (2.29). By using the result (3.45) for a, we arrive at the following 
expression for cr/S: 



S 2'kT 



11 9 2 n 

1 - - 2ci + I6C2 + -It 
6 3 



115 

— Ci - 32c2 + 12(2c3 + C4) 



+ 0(/2^) . (3.48) 



It was conjectured that e^ /{2'kT) would be the lower bound for a jr. in [20], at least when 
/i = 0. However, this possibility must again be questioned in light of the fact that the universal 
behaviour observed to leading order is again affected by the four-derivative couplings. 
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Finally, let us turn to the thermal conductivity (3.29). The following interesting ratio 
was constructed in [17]:^^ 

= 47r2 (^1 + ^(23ci + 24c2) + ^(37ci + 36(c2 + Acs + 2c,))lf + 0(/i^)^ . (3.49) 

Here we note that again the four-derivative interactions modify the leading behaviour but, in 
particular, also introduce a dependence on jl, similar to what was found for the ratio r//s in 
(3.25). 

4. Discussion 

In this paper, we calculated the thermal and hydrodynamic properties of the CFT plasma dual 
to a charged planar AdS black hole. These calculations were made within the perturbative 
framework where the leading Einstein-Maxwell action was extended to include a general set 
of four-derivative interactions (2.1). We should say that this analysis partially overlaps with 
previous results and so let us briefly summarize what was already known in the literature: 
In [40], the authors considered the transport properties of the charged planar AdS black 
hole solution with the standard two-derivative action, including the electrical conductivity a. 
Our results agree with theirs. Born-Infeld black holes were considered in [41], i.e., the two- 
derivative action was extended to include the combination of four-F terms arising from the 
expansion of the DBI action. They found that r//s = l/47r which is obtained straightforwardly 
from (3.25) by observing that the ratio is independent of Cg and C4. The eflcect of adding a 
Gauss-Bonnet term to the gravitational action was considered in [42] and -q/s was calculated 
in the charged planar AdS black hole background. In [43], this was extended to include 
both the Gauss-Bonnet term and the two independent four-_F terms. In both cases, our 
results are in agreement with theirs. The work was also extended in [44] to include a dilaton 
coupling to the Gauss-Bonnet term. However, in as discussed in [14], this extra coupling 
does not effect our results with the present perturbative approach. The authors of [45] had 
previously considered the thermodynamics with Gauss-Bonnet gravity and the iom-F terms 
for charged AdS black holes with flat, spherical and hyperbolic horizons. Our results for the 
thermodynamic behaviour agrees with theirs for the flat case. Finally in [46], the effect of 
the RabcdF"-'' F"'^ interaction was considered on a when = 0. Their result is reproduced 
by setting Ci = = fi in (3.45). Hence our comprehensive analysis agrees with all previous 
results where it should. 

As noted in subsection 2.3, our result for the entropy density (2.37) agrees with the 
result found using Wald's formula for higher curvature theories [36]. As expected then, the 
'geometric' expression for the entropy only involves Ci and C2 since it is only for these couplings 
that the corresponding four-derivative interactions in (2.1) involve the curvature tensor. It 

^■^Note that the leading order result given in [17] was Stt^ because their normalization for the gauge kinetic 
term differs by a factor of 2 from that used here. 
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is a nontrivial result that the couphngs C3 and C4 still do not enter when this result is re- 
expressed in terms of the temperature and chemical potential. We might note that these 
two coefficients appear in (2.33) and (2.35) and so a nontrivial cancelation is required for 
the final result in (2.37) to be independent of these parameters. It is interesting to observe 
that a similar cancelation occurs for C2 when the entropy is expressed in terms of T and /x, 
so that s is also independent of this parameter. Referring back to (2.4), this means that the 
entropy density is the sum of two contributions proportional to each of the central charges 
a, c appearing in the CFT 



where g{n = 0) = 1 = h{iJ, = 0). The shear viscosity also depends only on Ci and C2, as 
can be be inferred from (3.24) or (3.25). In this still appears in the result when it is 

expressed in terms of T and /x but again Cg and C4 do not appear. That is, the shear viscosity 
in the CFT depends on the central charges a, c but also the coupling of the stress tensor to 
the U{1) current parameterized by C2. 

As noted above, only Ci and C2 appear in corrections to the shear viscosity. This result is 
in keeping with the spirit of the recent work in [47]. This work considers generalized higher 
curvature theories of gravity and attempts to relate the shear viscosity to a 'gravitational 
coupling' evaluated at the horizon with a Wald-like formula [36]. While the conjectured 
formulae reproduce known results for certain higher curvature actions, it is known that these 
expressions fail to reproduce the correct shear viscosity in complete generality [38, 48]. 

All four of the couplings, Ci, C2, and C4, appear in our expressions for the charge density 
(2.38) and conductivity (3.45), as well as the free energy and energy densities. Since C3 and C4 
parameterize couplings in the four-point function of the U{1) currents, it is natural that they 
play a role in correcting the properties of the CFT plasma directly related to the corresponding 
charge. It is interesting to note that these two couplings only appear in (2.38) and (3.45) in 
the combination 2c;, + C4. In fact, examining all of our expressions in sections 2 and 3, one 
finds that it is only this particular combination of C3 and C4 that appears everywhere. Hence if 
we organized the four-F interactions in the action in terms of (F^)^ — 2F* and, say, F'^, then 
the first term would completely decouple from the present analysis. This observation was 
previously noted in [45] . We should add that this combination seems only to be distinguished 
by the particular black holes that we are considering here. For example, this combination 
does not appear as the four-F interaction in the low-energy expansion of the Dirac-Born-Infeld 
action [49] or in the four-derivative extension of the N = 2 supergravity action [26] . 

In the extremal limit (2.34), the entropy density (2.37) becomes 



s = —cT^ g{li) - — 



aT^h{fi), 



(4.1) 



s = 




(l-48(ci + C2)) , 



ii^ (1 - ci - 24c2) 



(4.2) 



3\/64^ 



- 22 - 



This result reflects the fact that the extremal black hole still has a finite size horizon. The 
interpretation in terms of the dual gauge theory is that even at zero temperature, a finite 
chemical potential will produce a deconfined 'plasma' in the dual CFT. It is interesting to 
note that the ratio r]/s for this extremal plasma, given in (3.26), only depends on C2. We 
should say that we expect that this 'exotic' behaviour of the CFT at zero temperature reflects 
our restriction of including only the metric and a single vector in the gravitational theory. 
That is, in many supergravity scenarios, the gauge kinetic terms will couple to various scalars 
and from the dual CFT perspective, such couplings reflect a nontrivial three-point function 
mixing two currents with some scalar operator. In such a scenario, the area of the horizon 
typically shrinks to zero size in the extremal limit [50] and so a finite temperature would be 
required to produce a deconfined plasma. 

In our general analysis, the coefficients q are treated as independent couplings. As 
described above, N = 2 supergravity in five dimensions provides a particularly interesting 
class of theories, since the super-graviton multiplet also contains a U{1) vector. In this case, 
the bosonic action for the metric and this vector takes precisely the form given in (2.1) 
with K = l/A^/S. However, in this case, the bulk supersymmetry is sufficiently restrictive to 
constrain all of these four-derivative couplings to be proportional to a single overall constant. 
Given the result in (2.4), this means that all of these coefficients are proportional to (c — a)/c 
where a and c are the central charges of the dual super symmetric CFT. In this case, the 
vector in the supergravity muliplet is dual to the CFT's i?-current, which is in the same 
supermultiplet as the stress tensor [27]. Hence the previous result is in agreement that the 
observation that two- and three-point functions of these two operators in the CFT (at zero 
temperature and chemical potential) are parameterized entirely by the two central charges 
[51]. 

After examining the supergravity action [26] in more detail in appendix A, we find that 

1 1 c — a , , „, 

■== = -^■='-16— ■ 

Hence for these theories, from (3.25), the ratio ij/s becomes 



rj 1 

s A-K 



1-8C,- ^2^'^^ 



3(1 + ^1 + 2/273)' 



(4.4) 



in the presence of a chemical potential. It is clear that the sign of the third term is con- 
trolled by Cj and in fact, this sign will be the same as that appearing in the second term. 
Hence if Ci is positive, both of these contributions lead to a violation of the conjectured 
KSS bound [16]. So in this particular class of theories, introducing a chemical potential 
only makes the violation stronger. For example, we note that for large Jl (4.4) yields: 
•q/s l/(47r) (1 — 24ci -|- 0(ci//i^)) . We should also add that in fact it was found that c> a 



^^A similar result appears in [53]. 
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and hence Ci > for all of the examples of superconformal gauge theories examined in [14]. 
Hence it appears that such violations of the KSS bound should be considered generic rather 
the exception to rule. We return to this point below. 

It is also interesting to examine the bounds conjectured in [20] for the conductivity in 
the case of supcrgravity. If wc restrict our attention to fx = 0, as was explicitly considered in 
[20], (3.46) and (3.48) reduce to 



rj 



10 , , 

= l-—c, + IQc, , (4.5) 

/i=0 

= -|- (1 - 2c, + 16c,) . (4.6) 

Hence while the general expressions also depended on and C4, only Ci and c, appear in 
both of these ratios when jl vanishes. Now let us consider these results when we substitute 
the supergravity result (4.3), c, = — Ci/2, and assume that is positive, as found in specific 
examples [13, 14] but is plausibly a general result. In this case, the higher order corrections 
reduce the value of both of the ratios, (4.5) and (4.6). Hence the first result is in agreement 
with the conjecture that (TT^/(ry e^)|^^Q < 1. However, the second ratio was conjectured to 
obey a lower bound cr/S |^^q > e^/{2iTT) and so the correction produces a violation of this 
conjectured bound. 

The fact that the sign of contribution coming from the chemical potential in (4.4) was 
controlled by is related, in part, to the fact that only even powers of the chemical potential 
everywhere in our analysis — of course, Uq has an overall factor of This property arises 
naturally from the tensor structure of gravitational action (2.1) and the particular background 
that we are studying, i.e., all of the relevant interactions contain even powers of the field 
strength Faj,. Of course, the bulk vector appears with an odd power in the two Chern-Simons- 
like terms in (2.1) but these interactions did not play a role in the present calculations, e.g., 
the two couplings, k and Cg, appear nowhere in our results. However, the two derivative 
coupling K is known to play a role when rotation is introduced [52]. Further, these couplings 
should play a role if one considers the magneto-hydrodynamics of the CFT plasma, i.e., if 
we introduce both bulk electric and magnetic fields, as has been recently studied with the 
AdS/CFT correspondence for three-dimensional field theories [54]. Extending this work to 
four-dimcnsional CFT's as considered here may be an interesting direction for future research. 

On the other hand, the behaviour noted above suggests that the properties of sQGP 
studied at RHIC, with (n^/T)'^ ^ 0.02, should be almost unaffected by the baryon chemical 
potential. Of course, one may question whether or not this property of our holographic 
models carries over to the sQGP. However, it seems that the appearance of only even powers 
of fi should be a general feature emerging from the CPT invariance of the underlying gauge 
theory, irrespective of whether the latter has a holographic dual or is even conformal. For 
example, the behaviour of the plasma with a net 'quark' density must be identical to that of 
the charge-conjugate plasma with a net density of 'anti-quarks'. Hence, one should expect 
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that the thermal and hydrodynamic properties of the gauge theory plasma should generally 
be even functions of /j.. 

There was a striking difference in the computations in subsections 3.1 and 3.2. In par- 
ticular, the shear viscosity could be framed in terms of quantities which were independent of 
the radius and so the latter could be evaluated at the horizon. Hence the shear viscosity of 
the CFT could also be interpreted as the shear viscosity associated to the stretched horizon 
in the membrane paradigm [55]. In contrast, the quantities determining conductivity evolved 
nontrivially in the radial direction and so the same connection could not be made to the 
membrane paradigm. 

Our analysis of the conductivity contrasts with the discussion in [6] which considered the 
conductivity with a vanishing chemical potential. Note that in this case, the mixing observed 
(3.31) vanishes and one may set /it^ = 0. Further with /x = 0, as seen in (3.33), the effective 
mass M{u) vanishes in the low-frequency limit and so the radial evolution (3.35) becomes 
trivial. Hence, in the absence of a chemical potential, there is a simple relation between the 
conductivity in the CFT and the universal conductivity of the stretched horizon, a^^, = e^/g^ 



Since the conductivity is a dimensionful quantity, the factor a/ 9zz{uo) appears to convert the 
length scale in the CFT to the corresponding proper length at the horizon. Of course, our 
conductivity (3.44) simplifies to reproduce this result in the limit that q (or /j.) vanishes. 

It is interesting to consider the radial flow found with 7^ by 'evaluating the conduc- 
tivity' at an arbitrary radius, i.e., removing the limit u — from (3.39). This yields 



where in the second expression, we have evaluated the result for the leading order background 
with (3.43) and denoted the conductivity in (3.44) Hence the radial evolution is such 

that cr{u) decreases monotonically as the radius varies from n = 1 at the horizon to u = at 
the asymptotic boundary. We also note that the boundary condition at horizon is precisely 
a{u = 1) = crcFT,;i=o = Cmb fo/L. Hence the membrane paradigm still sets the inner boundary 
condition for the nontrivial radial evolution but in general then, CTcft < o'mhfo/L, where the 
equality is only achieved when /x = and there is no radial evolution. 

In general, our discussion in subsection 3.2 generalizes the discussion of [6] to include 
both a finite chemical potential and the effect of higher derivative interactions. While the 
simplicity of the shear viscosity computation presented there is essentially not effected by 
these additional complications, there are in fact additional simplications of the computation 
beyond our presentation in subsection 3.1. We will examine these further in an upcoming 
paper [48], as well giving a covariant Wald-type formula for rj. 

To close, we would like to comment on a possible implication of the 'gravity as the 
weakest force' conjecture, i.e., the recent conjecture in [21] that there should be a general 



[6]: 



O'CFT,/i=0 — CTjnb 



(4.7) 




(4.8) 
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upper bound on the strength of gravity relative to gauge forces in quantum gravity. This 
conjecture requires that there are always light 'elementary particles' with a mass-to-charge 
ratio smaller than the corresponding ratio for macroscopic extremal black holes and so allow 
the extremal black holes to decay. Recently, there has been some discussion of this conjecture 
in the context of the AdS/CFT correspondence and the implications for the spectrum of the 
CFT [56, 57] — see also footnote 4 in [20]. The relation which we would like to draw relies 
on the further corollary that higher derivative corrections should reduce the mass-to-charge 
ratio of extremal black holes in a consistent theory of the quantum gravity [21, 58]. In the 
present context, this suggests computing the ratio of the energy density to charge density 
PE/nq in the extremal limit (2.34) 

^ = f ^1 U - l!^i±i^^ . (4.9) 

(\ 3\/3 T 
— ) = — i=~To, is the leading order 'classical' result. Hence the weak 

gravity conjecture would impose the constraint 47ci +48c2 > 0. However, if we again consider 
the specific case of supergravity with C2 = this constraint becomes simply Cj > 0. 

While the above analysis is suggestive, it misses the intent of the discussion in [21, 58] 
which was phrased in terms of extremal black holes in asymptotically flat space. Their natural 
assumption was that the 'quantum' corrections coming from higher derivative terms in the 
gravitational action should decrease the mass-to-charge ratio but also do so in a way that 
the decrease becomes more pronounced for smaller (extremal) black holes. Hence a proper 
comparison requires repeating the analysis of section 2.2 for charged AdS black holes with 
spherical horizons. This is a straightforward exercise and the final result replacing (4.9) is 



EE = ( EE] ( 



7 + 34-^ + 107 + 138^ 
l-c,/(ro/L)) with/(ro/L)= ^ (4.10) 



Above To is the position of the horizon in coordinates where the area of the spherical horizon 
is 27r^rQ. In fact, a full analysis yields a ratio which depends on all of the Ci but in presenting 
(4.10), we have focussed on the supergravity case where all of these dimensionless parameters 
are proportional to Ci. One easily verifies that in the limit of large black holes, i.e., Tq/L S> 1, 
this result reduces to that for the planar black holes given in (4.9) when C2 = — Ci/2. It is 
also evident that /(ro/L) is positive for all values of r„/L and hence the higher derivative 
corrections always reduce the mass-to-charge ratio of these extremal black holes as long as 
Ci > 0. As desired, this effect is also largest for small black holes because of the factor 
of rljl? in the denominator of /(ro/i^). However, it is interesting that ji^^jV) does not 
decrease monotonically as ro/L grows. Instead the function exhibits a local minimum near 
To ~ L, which seems to be an effect of the asymptotic AdS geometry. 

In certain cases, it is well understood that there are supersymmetric bound states of 
giant gravitons [59] carrying the same charges as these charged black holes with spherical 
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horizons. Further that there is a gap in the spectrum between the extremal black hole and 
these bound states [60] . Hence one can understand the details of realizing the 'gravity as the 
weakest force' conjecture within this framework. However, our interest in these issues comes 
rather from the general constraint Ci > 0. In particular with (2.4), we can infer that the 
weak gravity conjecture requires an inequality for the central charges of any four- dimensional 
conformal field theory with a gravitational dual, namely c > a. Of course, this is precisely the 
inequality which was observed for the broad class of superconformal gauge theories examined 
in [14]. Hence it seems there may be a deep connection between this provisional observation 
for superconformal gauge theories and the consistency of their holographic duals as a theories 
of quantum gravity. 
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A. Field Redefinitions and the four- derivative action 

In this appendix, we will demonstrate that our action (2.1) contains the most general four- 
derivative interactions involving a single Maxwell field. We should add that a similar analysis 
appeared in [30] but they began with a more restricted starting point. First, we consider the 
leading order two-derivative action 



Next the most general four-derivative action for gravity coupled to a single U{1) vector takes 
the form: 



+ p^R'^'^FacFi,^ + PsRabcdF^'^F'"^ + PiRabcdF'^F'^'^ 

+55VaFbcV^F'^ + SeV^Fab F"^ + •JrVaV^'Ffc, F"^ + S^V'^VaFbc F"" 

where, as in the main text, we use = FabF"''^ and = F"-bF^cF'^dF'''a- Here all of the 
coefficients, cti, Pi, 5i,7i and K2 are dimensionless constants, that we expect are generically 




(A.1) 




(A.2) 
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very small. Many of the four-derivative terms above can be eliminated by simply integrating 
by parts. For example, 



/ 



C Z7 b 
C 



(A.3) 



= j d^X^ (VaF^^V^F"^ - R'^FacFi,'' + RabcdP'^'F^'^) ■ 



Using integration by parts, as well as the identities V^aFbc] = = R[abc]di o^i^ can eliminate 
(^4,5,6,7,8 and 72,3. In this way, the general four-derivative action can be reduced to 



d^Xy 



aiR^ + a2RabR"'^ + Oi^RabcdR"^'^'^ 



(A.4) 



+ p^R^'^FacF,,^ + PsRabcdF^'^F^'^ + 5i {F^f + 52F^ 

+5^VFabV''Fj> + Jie'^^'^'FabF^d'^fFfe + K2e''^^'AaRbcfgRde^' 

Now consider making field redefinitions: gab — gab + ^9ab and Aa ^ Aa + SAa- The most 
general field redefinition involving two-derivative contributions can be written 



Sgab = HiL'^ Rab + l^2L^ FacFb" + {n^L^R + haL'^F^ + //g) gab , 

SAa = \lAa + \2^^Fba + abode F^'F'^^ . 



(A.5) 



Note that and Ai give a (constant) rescalings of the metric and vector, respectively, which 
will prove useful in the following. We also note that other than the rescaling the field redefini- 
tion of the vector involves covariant terms so that the modified action remains invariant with 
standard gauge transformations. In general, we might note that, aside from the two rescal- 
ings, the field redefinitions (A.5) contain six two-derivative terms while the four-derivative 
action contains eleven interactions. Hence, on general grounds, we expect that wc will be left 
with five independent terms in our higher-order action. With the field redefinitions (A.5), the 
leading change in the action comes from the variation of (A.l) 



5L 



1 



2£p 



6 



1 



1 



9 



,ab 



_|_ ^jpcicpib 



Sg, 



'ab 



+ 



(VfeF^" + Kie'''""^'FbcFde) M„| . (A.6) 



Of course, one should note that we have integrated by parts to produce the expressions in 
(A.6). Now we will divide this variation into two parts, examining separately the contributions 
to the four- and two-derivative actions. Beginning with the former, one finds: 

^ ^^{^ + lns)R'-^^iRabR''' + 



SL 



d^Xy 



-9 



1^1 



f - f + ) RF-(A.7) 



+ 



(f - ^2) R^'Fa^F,^ + - Y + A3) {F-f 



+ [y- ^^''^^3} + A2 V^F^bV^i^c' + {1^1X2 + A3) e^'''''^FabF''''V^Ffe 
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An obvious choice of the field redefinition parameters is then 



Ml = "2, tJ'2 = (32 + (X2/2, fj,3 = -{2ai + q;2)/3, 

/X4 = - (24/3i - 12/32 + 2ai - llaa) /36, A2 = -S3, A3 = -71 + K1S3. (A.8) 
These choices ehminate six of the four-derivative interactions leaving 

L2 



I d^x^ [ a^RabcdR"''"^ + PsRabcdF^'^F^ (A.9) 



where 



^1 = (^1 + ^ (2ai - 29a2) + 3^ (/^i - /?2) - 8^171 + SnlSs , 

^2 = ^2 + Y + y + 16^^171 - 16k? <53 . (A.IO) 

One may try to be more clever with the field redefinitions in the present case. In particu- 
lar, as seen in [45] (or our calculations in sections 2 and 3), when the background only contains 
a radial electric field, the four-F terms only couple to the graviton and gauge field equations 
for motion with the combination 2Si + 62- Hence might try to make a field redefinition that 
sets this combination to zero. If we go back to (A. 7), we find that the field redefinitions yield 
the following variation of this linear combination: 

A [2Si + S2] = \ {^12 - l^i) . (A.ll) 

Hence we can indeed arrange to set this combination of couplings to zero. However, the result 
will be that wc would not be able to eliminate all of the R^ and RF^ interactions involving 
the Ricci tensor and Ricci scalar. Hence such a field redefinition will simply replace the 
complications of accounting for the four-F interactions with those of accounting for another 
set of interactions. With respect to the dual CFT, these two four-F interactions define two 
characteristic parameters that would appear in the four-point function of the dual current. 
So it seems to natural not to field redefine them away, as the four-point function must be 
invariant and the previous parameters would just appear from exchanges between bulk three- 
point interactions rather than being manifest in a four-point contact interaction. However, 
given this discussion, we could arrange the four-F terms in (A.9) as 



^3 / d'-V^9 [I (2^1 + ^-2) {F^r + ^2 {f^ - \ [F-f) 



(A.12) 



in which case, the second term should not contribute in the present calculations. 

At this point, we have not yet commented on fixing the values of ^5 and Ai. In this 
regard, let us consider the effect of the field redefinitions on varying the two-derivative action 
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(A.l): 



5U 



12 5 / 3 



(A.13) 



^ _ A 
8 2 

Hence the convenient choice which we make in fixing these parameters is to set 



/X5 = -4 {ill + 5//3) , Ai = //I + 12//2 + 5//3 + 60//4 . 



(A.14) 



This choice leaves the Planck scale fixed, as well as the coefficient for the vector kinetic term. 
Hence after the field redefinitions, the two-derivative action (A.l) becomes 



1 



lip 



12 



de 



(A.15) 



where 



-1 = J_ (1-10 (^1 + 5^3)) 

Kl = Kl (1 + 3 {fii + 12/Lt2 + 5/X3 + 6O/X4)) . 



(A.16) 



Next we would like to apply this analysis to compare the four-derivative supergravity 
action given in [26] to the effective action (2.1) studied in the present paper. Of course, 
the leading two-derivative terms in the supergravity action take precisely the form given in 
(A.l) with Kl = 1/4^3(1 — 32ci) - here and in the following, we present the results in terms 
of Ci = (c — a) /{8c), as in (2.4). The four-derivative supergravity action can certainly be 
described in terms of our general action (A. 2) where the dimensionless coefficients are all 
assigned specific values proportional to Ci. However, the supergravity action is naturally 
described in terms of the Weyl tensor and so to make this matching, we must first re-express 
the following: 

f^abcd T) ryabcd ^ jj nab i ■'■ r?2 /a i 7\ 

'-'abcd^ — ^abcd^ — -^^ab^ 6 ' \-^-^') 

CabcdF'^F'"^ = RabcdF'^F'"^ - ^R'^FacFb'' + IrF^ , 

3 6 

which applies for the five-dimensional Weyl tensor. Now as described above, integration by 
parts can be used to eliminate many terms in the general action (A. 2), reducing it to the 

^^Note that this comparison requires care since [26] adopts the supergravity conventions of [28], with, e.g., 
the mostly minus convention for the signature of the metric. We thank Sera Cremonini for her detailed 
explanation of these conventions. 
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form given in (A. 4). Putting the supergravity action in this particular form yields: 



1 




4 




02 = 




1 


(32 = 


4 






41 




5 


"288'^'' 


52 = 




5 




1 




K2 = 





as = Cl , 

/9i = 7Ci, /?2 = -^Ci, /33 = -lci, (A.18) 

<^3 = 2Ci , 

7, = ^ci , n2 = ^ci . (A.19) 

Now applying the field redefinitions (A.5) with the parameters fixed as in (A.8), we are able 
to further reduce the four-derivative action to the canonical form (A. 9). While these field 
redefinitions leave the values of 03, and K2 unchanged from those given above, (A. 10) 
yields 

5i = ^ci, ^2 = -^Ci. (A.20) 

As a final step, we give the Einstein and Maxwell kinetic terms their standard normalization 
with (A. 14) yielding 

1 1 10 \ . , 1 , , , 

L^^I?\ ' '^i = ^1 (1 - 256ci) = ^ (1 - 288ci) . (A.21) 

Hence, after all of these manipulations, the supergravity action takes the form given in (2.1) 
with 

_ 1 _ 1 _ 5 _ 1 

C2 — 2*^1 ' — 2^Ci , C4 — — Cl , Cg — ^^^^1 ' (A. 22) 

as well as «; = 1/4^3(1 - 288ci) and Ci = (c-a)/(8c), as in (2.4). As discussed, this latter 
combination of the central charges in the (supersymmetric) CFT is seen to explicitly fix the 
coefficients all of these four-derivative corrections to the supergravity action. 
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